Nematic phase and phase separation near saturation field in frustrated ferromagnets 
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We argue the effect of quantum fluctuations on the magnetization process of the quantum frus- 
trated ferromagnets. It is found that on general grounds in the neighborhood of the ferromag- 
net/antiferromagnet classical lst-order phase boundary in zero external field, the phase separation 
or the non-classical phase must appear slightly below the saturation field in the quantum case, if 
the classical AF is not the eigenstate. Then, we study the ferromagnetic J\-J<i S = 1/2 Heisenberg 
model (Ji < 0) on the bcc lattice from the viewpoint of the magnon Bose-Einstein condensation. 
For -1.50097 < Ji/J 2 < -1.389, the nematic phase is expected and for -1.389 < Ji/J 2 < -0.48 
the phase separation appears under high magnetic field. 
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Quantum fluctuations (QF) in frustrated magnets can 
introduce a drastic change in the ground state from the 
classical system and induce various exotic phases. As one 
of quantum phases, the spin nematic has been recently 
paid much attention, which is nonmagnetic but exhibits 
the long-range order of the multiple spins. 

In this paper we discuss the effect of QF on frustrated 
magnets which have dominant ferromagnetic (FM) cou- 
plings, studying the possibility of the spin nematic phase. 
We consider the S = 1/2 frustrated Heisenberg model de- 
scribed by 



H 
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The classical limit (S — »• oo) is equivalent to considering 
the spin operator as the vector of the fixed length S. In 
the case with one magnetic ion per unit cell, the classi- 
cal ground state is obtained by minimizing the Fourier 
transform of the exchange interactions 
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The Ji- J 2 model, described by the nearest-neighbor cou- 
pling J\ and next-nearest-neighbor coupling J2, has been 
extensively studied in the research of frustration effects. 
In the one-dimensional (ID) J1-J2 chain with FM J\ < 0, 
the spiral spin structure is classically stable for | J\\/ J2 < 
4 and, is replaced by the spin nematic phase in S = 1/2 
system in high magnetic field 1 . In LiCuVCU, the mag- 
netic properties of which is understood as the quasi- ID 
J\- J 2 model with the weak interchain couplings^, a non- 
trivial phase is recently found slightly below the satu- 
ration field and it is expected that the spin nematic is 
observed for the first time 3 . 

The J\- J 2 model on the square lattice is another sys- 
tem which exhibits the nematic phased. Classically, three 
types of orders are stabilized according to J\^'- Neel anti- 
ferromagnetic (NAF), collinear antiferromagnetic (CAF), 
and FM phases. In the quantum case, the nematic phase 
was proposed from both the numerical and analytical 



approaches^ near the classical FM/CAF phase bound- 
ary (J1/J2 = —2 and J\ < 0), though there is still de- 
bate about the stability of spin nematic order in the case 
of zero field££ Slightly below the saturation field, the 
spin nematic phase is firmly induced by the two-magnon 
instability in the saturated FM state for the broad pa- 
rameter range —2.5 < J1/J2 ^ — 0.2i£. The recent ex- 
perimental advance found the various compounds, e.g., 
BaCdVO(P04)2, which are adequate to consider as the 
FM J1-J2 model-. Some of them lie in the parameter 
range where the nematic phase is theoretically proposed. 
It would be useful to discuss how QF changes magnetic 
properties near the boundary between FM and AF phases 
in general frustrated FM systems. 

In this paper, we first argue the quantum effects near 
the FM/AF phase boundary on general grounds. In the 
neighborhood of the FM/AF- 'classical' phase boundary 
in zero field, there necessarily exists either a non-classical 
phase or a phase separation accompanied with magneti- 
zation jump below the saturation field, if the FM/AF 
transition is first order and the classical AF state is 
not the exact ground state of the Hamiltonian. As the 
FM/CAF phase boundary in the S = 1/2 square- lattice 
J\- J 2 model satisfies this condition, this perspective rein- 
forces the former proposal of the quantum spin nematic 
phase. As another example of this argument, we further 
study the three-dimensional J\- J 2 model on the bcc lat- 
tice, which has been extensively studied as one of the 
simplest frustrated magnets. Recent theoretical studies 
concluded that in the AF case (Ji,2 > 0) the classically 
expected AF orders persist even in the quantum case£ 
For the FM J\ case, however, the above discussion as- 
serts the appearance of non-classical behavior. In the 
later half of this paper, we study the appearance of spin 
nematic phase or the phase separation in high field in the 
S = 1/2 J1-J2 model with FM J 1 on the bcc lattice. We 
employ the dilute Bose gas approach from the viewpoint 
of the magnon Bose-Einstein condensation (BEC) 9 . 

Let us start with a general argument on the quan- 
tum effect on the phase boundary between FM and AF 
phases. In the absence of magnetic field, the classical 
ground state usually shows a first-order phase transi- 
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FIG. 1: (color online) Three-dimensional body-centered cu- 
bic (bcc) lattices. Filled circles denote spins connected by 
Heisenberg exchange interactions. The bcc lattice is bipartite 
when J2 = and can be divided into two sublattices, which 
are distinguished by the size of spheres. 



tion between FM and AF phases. In most settings of 
the quantum systems, the classical AF state |AF;CL) 
is not an eigenstate of the Heisenberg Hamiltonian, 
#|AF;CL) = £ C l|AF;CL) + | a ), where (AF;CL|a) = 
and (a\a) > 0, albeit the FM state |FM;CL) is the eigen- 
state having the same energy as the classical one. The 
variational principle guarantees that the true quantum 
(qAF) ground state in the classical AF regime has the 
lower energy than that of the FM and AF Neel states 
at the classical FM/AF phase boundary. Therefore, 
FM/qAF phase boundary shifts from the classical one 
toward the FM phase by a certain amount in zero field. 

The magnetization process in applied field must be 
compatible with this fact. However, if one considers the 
saturation field given by one-magnon flips, which is al- 
ways the same as the classical value H c , the saturation 
field must vanish at the classical FM / AF phase boundary. 
This concludes that the true saturation field is not given 
by single magnon instability. The AF phase described 
with the single magnon BEC below the saturation field 
is veiled by qAF phase near the FM/AF phase bound- 
ary. We thus conclude the appearance of a non-classical 
phase or the phase separation below the saturation field. 
An appealing alternative to the single magnon BEC in 
applied field is the BEC of multiple-bound magnons. 

Hence, as in the square-lattice J\-J 2 model, the QF 
plays a considerable role in the magnetization process in 
J\- J 2 models on other lattices, even in three dimensions 
where the QF is believed to be weak. The J\-J 2 model 
on the bcc lattice satisfies the condition discussed above. 
We hereafter discuss the possibility of the nematic phase 
and the phase separation in the magnetization process 
on the bcc lattice. 

Before studying the quantum effect, let us briefly re- 
view the ground state in the classical case. When the 
external field H is absent, the three types of phases ap- 
pear; NAF phase for J ± /J 2 > 3/2 and J x > 0, CAF 
phase for -3/2 < J 1 /J 2 < 3/2 and J 2 > 0, and FM 
phase for J1/J2 < —3/2 and J\ < 0, as shown in Fig. [2j 
In CAF phase, the lattice is divided into two sublattices, 
and the spins form the NAF on each sublattice, where 
the spin angle between two sublattices is independent 
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FIG. 2: (color online) Spin configurations for the three phases 
('FM', 'NAF' and 'CAF') in the bcc lattice, (i) 'FM' repre- 
sents the fully polarized phase, (ii) In 'NAF', the spins on 
each sublattice align ferromagnetically while those on differ- 
ent sublattices are anti-parallel, (iii) 'CAF' is made up of two 
antiferromagnetically-ordered sublattices, which, as a whole, 
align in a collinear manner. 



and QF fixes it. When the external field is applied in 
the CAF phase, the spins gradually point upward from 
the plane perpendicular to the external field. The angle 
between spins on two sublattices is independent even in 
the magnetization process. Within the linear spin wave 
theory, the (deformed) CAF phase is stabilized 8 . 

We can apply the argument discussed above to this 
model. At the FM/CAF phase boundary J 1 /J 2 = -3/2, 
the classical CAF state is not an eigenstate in the quan- 
tum Heisenberg model. Hence, near J\/J 2 = —3/2, a 
nontrivial quantum phase or the phase separation must 
appear under high magnetic field. On the other hand, 
at the FM/NAF phase boundary Ji = 0, the NAF state 
is an eigenstate even in the quantum case and then a 
quantum fluctuation is not important. We hence focus 
our attention on the CAF and FM phase boundary. In 
particular, we study the magnetic behavior slightly be- 
low the saturation field in a fully quantum manner using 
the method of a dilute Bose gas^. 

The emergent magnetic order below the satura- 
tion field can be viewed as a condensation of the 
magnons, which successfully explains various experimen- 
tal results*^. In the hardcore boson language, the spin 
operators at site I are written as^ Sf = —1/2 



a\ai, 



— a j\ an d Sf = a/, where a% represents the anni- 
hilation operator of a hardcore boson. The vacuum |fi) 
corresponding to the saturated FM state is defined by 
ai\Q) = 0. The Hamiltonian reads 



H = ^(cj(q) - /i)a^a q - 
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with a;(q) = e(q) - e min , /; = H c H. H ( . = e(0) - e min , 
and = 2(e(q) + U), where e(q) is given in Eq. (j2j), 
e m i n is the minimum of e(q = Q), and [/(—>• 00) is the 
hard-core potential. We see that the BEC occurs when 
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the external field is H < H c (fi > 0), which leads to 
(Sf) = (at) 7^ 0. The induced spin-ordered phase pos- 
sesses characteristic wave vector Q. To clearly specify the 
characteristics of BEC, we call the BEC induced by the 
magnons with the magnetic quantum number +1 as the 
'single' magnon BEC. In large spin S systems, where the 
theoretical method is extended by using the spin wave 
expansion^ the classical phase near the saturation field 
can be well described with the single magnon BEC. In 
the quantum case, the single magnon BEC can induce a 
different spin-ordered phase from the classical one, albeit 
the saturation field of the single magnon BEC is exactly 
the same as the classical one. 

Since the method is conventional and detailed in 
Refs. flolll3lll3 . we mainly discuss our results on the 
bcc lattice. The dispersion relation is given by e(q) = 
4Ji cos(q x /2) cos(q y /2) cos(q z /2) + J2(cosq x + cos^ + 
cosq z ) and takes minimums at ±Q where Q = (7r,7r,7r) 
for J1/J2 > —3/2. The thermal potential per site E/N 
of the dilute Bose gas is determined by the interaction 
among the bosons condensed at q = ±Q. Then, the den- 
sities of the condensed bosons p±Q in the ground state 
are obtained by minimizing E/N . If we introduce the 
renormalized interactions between the same bosons and 
that between different ones respectively as Y\ and Y 2l the 
energy density E/N in the dilute limit is given by 

§ = y (pq + p-q) + {r 2 + r 3 cos 2(0 Q - 0- Q )}p QP - Q 



p(pq + p-q), 



(4) 



where (a q ) = y/Np^ exp(i6 q ) and the term with T3 ap- 
pears from the lattice structure. If 1^ is positive (neg- 
ative), the phases 0±q are pinned as 0q — 6Lq = mir 
[(ra+ |)7r] with an integer m. T M are calculated from the 
spin Hamiltonian later. 

Various ground states can appear according to the 
values of Y^. If T 2 - |r 3 | > T 1 > 0, the ground 
state is given by pQ = p/Y\ and p-Q = (or vice 
versa), breaking the chiral symmetry. This phase is, 
however, not found in the present J1-J2 model in the 
following calculation. If Ti > Y 2 — 1 17*3 1 > — Ti, on 
the other hand, the two modes condense simultaneously 
with the same density and the ground state is deter- 
mined as pq = p_ Q = p = /i/(ri + T 2 - |r 3 |) and 

M = v / P{ ei(Q ' Rz+0Q) +e^- Q - R ^-Q)}. In this phase, 
the spin expectation values are given as 



4p COS 2 (Q-R; 
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For positive 1^3, the CAF appears, as expected within 
the linear spin wave theory. Meanwhile, for negative 1^3, 
the emergent phase has a modulation of the density (Sf ), 
which can be considered as the spin-supersolid phase. 

If Ti < or ri + Y 2 — \Y 3 \ < 0, the low-energy bosons 
around q = ±Q attract each other and we hence expect 



a first order transition, or equivalently, phase separation. 
So far, our argument based on Eq. dU) is restricted to the 
magnon BEC that occurs in the single-particle channel. 
However, strong attraction can also induce formation of 
magnon bound states, which lead to the BEC of magnon 
bound states. 

The condensed phase of the two-magnon bound states 
has the different properties from the single-magnon BEC. 
The striking one is the absence of the transverse lo- 
cal magnetization, and instead the existence of the long 
range order in the quadratic channels (S^Sj~) 7^ for a 
certain bond (i, j), which corresponds to the spin nematic 
order. 1 ? 4 The stability of the two magnon bound state is 
implied from the divergence of the scattering amplitude 
of the two magnons. 

To study r M concretely, we study the ladder diagram 
shown in Fig. [3j which describes the exact scattering am- 
plitude M on the saturated FM phase; 



M(A,K;p,p') = V-p + ^-p'- 



d 3 p" M(A, K; p, p")(V- p » + V- P >- P ») 



(2tt) 3 w(K/2 + p") + w(K/2 - p") + A 
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where A is the total energy and K is the center-of-mass 
momentum. are given by Ti = M(0, 2Q; 0, 0)/2, 
T 2 = M(0,0;Q,Q), T 3 = M(0,2Q;0,qi)/2 with qi = 
(0,0,27r). The dispersion relation of the bound state 
A(K) is identified as the divergence of M by scanning 
it with respect to A and K. By using the harmonic ex- 
pansion, Eq. ((6]) is exactly solvabl e) 10 ? 13 ! 14 . 

First, we discuss the result of obtained for A = to 
study the single magnon BEC. For -0.48 < J 1 /J 2 < 3/2, 
the CAF appears since T 1 > T 2 - \T 3 \ > -T 1 and T 3 > 
hold. For -1.389 < Ji/J 2 < -0.48, we obtain Ti < 
and the phase separation near the saturation field is 
expected. The resulting phase below the magnetization 
jump may be CAF, but we do not exclude the possibility 
of nontrivial phases such as the spin nematic appearing 
by the lst-order transition. At J1/J2 = —1.389, both Ti 
and T3 diverge, implying the appearance of bound states. 
For -3/2 < J x /J 2 < -1.389, we have T 1 > T 2 - \T 3 \ > 
—Ti and T 3 < 0, which suggests the spin supersolid. 

To see the possibility of the bound magnon BEC, 
we detail the binding energy of the bound state. The 
minimum of the dispersion A(K) is always obtained at 
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FIG. 3: Scattering amplitude M given by the ladder diagram. 
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FIG. 4: (color online) Binding energy A of magnon pairs with 
the center-of-mass momentum K = (27r, 2tt, 2tt). In case of 
positive A, the bound magnon BEC is expected. 
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pairs closes earlier than that of the single magnons and 
thus spin nematic ordering veils the expected spin super- 
solid phase. For -1.50097 < Ji/J 2 < -3/2, the ground 
state is not the classically expected FM phase in the ab- 
sence of the external field since the reference state (FM) 
is destabilized by the fluctuation of the bound-magnon 
pairs. Our finding shows that the nematic phase is the 
serious candidate of the ground state. From the symme- 
try perspective the FM/nematic phase transition in zero 
field may be lst-order as in the case of the square lattice^. 
We note that there remains the possibility of appearance 
of CAF phase in very low field regime accompanied with 
the magnetization jump below the saturation field since 
the energy of the bound state is almost the same as the 
single magnon. 

Finally, we discuss the properties of the bound magnon 
condensed phase. The wave function of the bound- 
magnon pairs with respect to the relative coordinate is 
given by 



x(p) ^ 



cos 
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FIG. 5: (color online) Critical external field H c . For 
H > H c , the spins form the saturated FM state. The 
straight line (blue) denotes the critical field where the gap 
of bound magnon pairs closes. The short dashed line (red) 
denotes the field where the single magnon gap closes. For 
-1.50097 < Ji/J 2 < -1.389, the bound magnon BEC is 
expected, 'lst-order' represents the phase separation. The 
vertical dashed line is a schematic phase boundary between 
the bound magnon BEC and the single magnon BEC. 



Ko = (27r, 27r, 2tt). Then, A is given by solving 
d 3 p 



Ji 



cos 2 2f cos 2 cos 2 ^ 



7T 3 ou(K /2 + p) + o;(K /2 - p) - A 



1 (7) 



withcj(K /2 + p)+cj(K /2-p) = 8J1 +2 J 2 (3 + cosp x + 
cospy + cosp z ) for J1/J2 < —3/2 and 6J2 — 2 ^(cosp^ + 
cos p y + cosp z ) for J1/J2 > —3/2. The total bound- 
magnon energy is given by Eb = —2fi — A. 

The estimated results of A and the critical external 
field H c are respectively shown in Fig. [4] and Fig. [5] For 
—3/2 < J1/J2 < —1.389, the gap of the bound-magnon 



where we have omitted the normalization constant. The 
nematic order parameter, which is dominant on the 
nearest neighbor bonds, has the /a^ 2 -wave symmetry, 
i.e., (S- # i±(ex+ey+ez)/2 ) = 



), with the 



(StSi±(e x -e y -e z )/2) 

wave vector Ko. 

In conclusion, we have discussed the effect of QF on the 
quantum frustrated ferromagnets under high magnetic 
field. In general, near the FM/AF classical lst-order 
phase boundary, the non-classical phase or the phase sep- 
aration must occur in magnetization process, if the clas- 
sical AF is not an eigenstate. We stress that the phase 
separation occurs even in the isotropic Heisenberg model 
and is purely the quantum phenomena. As a concrete 
model satisfying the above condition, we study the J±- J2 
model on the bec lattice using the method of a dilute 
Bose gas. We found that for -1.389 < Ji/J 2 < -0.48, 
the single magnon condensed phase is not stable near 
the saturation field and the phase separation occurs. For 
-1.50097 < J1/J2 < -1.389, the bound magnon BEC 
which leads to the nematic phase is expected. 

We thank K. Totsuka for discussions. 
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In Ref. Q, although discussed, the region of the nematic 
phase is not declared explicitly. And for —0.2 < J1/J2 < 0, 
we do not comment on the possibility of the nematic phase 
because of the problem of the numerical precision in our 
calculations. 



